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Abstract 

We present the GLoBES ("General Long Baseline Experiment Simulator") software 
package, which allows the simulation of long-baseline and reactor neutrino oscillation 
experiments. One part of the software is the abstract experiment definition language 
to define experiments with beam and full detector descriptions as accurate as possible. 
Many systematics options are provided, such as normalization and energy calibration 
errors, or the choice between spectral or total rate information. For the definition 
of experiments, a new transparent building block concept is introduced. In addition, 
an additional program provides the possibility to develop and test new experiment 
definitions quickly. Another part of GLoBES is the user's interface, which provides 
probability, rate, and A% 2 information for a given experiment or any combination 
of up to 32 experiments in C. Especially, the A% 2 functions allow a simulation with 
statistics only, systematics, correlations, and degeneracies. In particular, GLoBES can 
handle the full multi-parameter correlation among the oscillation parameters, external 
input, and matter density uncertainties. 



a Email: phuber@ph.tum.de 
b Email: lindner@ph.tum.de 
c Email: wwinter@ph.tum.dc 



1 Introduction 



Neutrino oscillations are now established as the leading flavor transition mechanism for 
neutrinos. In a long history of many experiments, see e.g. [1], two oscillation frequencies have 
been identified: The fast "atmospheric" and the slow "solar" oscillations, which are driven 
by the respective mass squared differences. In addition, there could be interference effects 
between these two oscillations, provided that the coupling given by the small mixing angle 
#13 is large enough. Such interference effects include, for example, leptonic CP violation. 
In order to test the unknown oscillation parameters, i.e., the mixing angle #13, the leptonic 
CP phase, and the neutrino mass hierarchy, new long-baseline and reactor experiments are 
proposed. These experiments send an artificial neutrino beam to a detector, or detect the 
neutrinos produced by a nuclear fission reactor. However, the presence of multiple solutions 
which are intrinsic to neutrino oscillation probabilities [2-5] affect these measurements. 
Thus optimization strategies are required which maximally exploit complementarity between 
experiments. Therefore, a modern, complete experiment simulation and analysis tool does 
not only need to have a highly accurate beam and detector simulation, but also powerful 
means to analyze correlations and degeneracies, especially for the combination of several 
experiments. In this paper, we present the GLoBES software package as such a tool. 

GLoBES ("General Long Baseline Experiment Simulator") is a flexible software tool to simu- 
late and analyze neutrino oscillation long-baseline and reactor experiments using a complete 
three-flavor description. On the one hand, it contains a comprehensive abstract experiment 
definition language (AEDL), which allows to describe most classes of long baseline and reac- 
tor experiments at an abstract level. On the other hand, it provides a C-library to process 
the experiment information in order to obtain oscillation probabilities, rate vectors, and 
Ax 2 - values (c/., Fig. HJ. In addition, it provides a binary program to test experiment def- 
initions very quickly, before they are used by the application software. Currently, GLoBES 
is available for GNU/Linux. Since the source code is included, the modifications to other 
operating systems should be doable. 

GLoBES allows to simulate experiments with stationary neutrino point sources, where each 
experiment is assumed to have only one neutrino source. Such experiments are neutrino 
beam experiments and reactor experiments. Geometrical effects of a source distribution, 
such as in the sun or the atmosphere, can not be described. In addition, sources with a 
physically significant time dependencies can not be studied, such as supernova?. 

On the experiment definition side, either built-in neutrino fluxes (e.g., neutrino factory) or 
arbitrary, user defined fluxes can be used. Similarly, arbitrary cross sections, energy depen- 
dent efficiencies, energy resolution functions as well as the considered oscillation channels, 
backgrounds, and many other properties can be specified. For the systematics, energy nor- 
malization and calibration errors can be simulated. Note that energy ranges and windows 
and bin widths can be (almost) arbitrarily chosen, including bins of different widths. To- 
gether with GLoBES comes a number of pre-defined experiments in order to demonstrate the 
capabilities of GLoBES and to provide prototypes for new experiments. In addition, they 
can be used to test new physics ideas with complete experiment simulations. Examples 
for these prototypes are the MINOS, ICARUS, and OPERA simulations from Ref. [6], the 
JHF-SK and NuMI superbeam simulations from Refs. [7,8], the JHF-HK superbeam upgrade 



1 



AEDL 
Abstract Experiment 
Definition Language 

Defines Experiments 
and modifies them 



| GLoBES | 




GLoBES User Interface 

C-library which loads 

AEDL-file(s) and 
provides functions to 
simulate experiment(s) 



Application software to compute 
high-level sensitivities, precision etc. 



Figure 1: General concept of the GLoBES package. 

simulation from Ref. [7], the neutrino factory simulations from Ref. [7], and the reactor 
experiment simulations from Ref. [9]. Other projects using earlier versions of the GLoBES 
software include Refs. [10-13]. 

With the C-library, one can extract the A% 2 for all defined oscillation channels for an 
experiment or any combination of experiments. Of course, also low-level information, such 
as oscillation probabilities or event rates, can be obtained. GLoBES includes the simulation 
of neutrino oscillations in matter with arbitrary matter density profiles, as well as it allows to 
simulate the matter density uncertainty. As one of the most advanced features of GLoBES, it 
provides the technology to project the A% 2 , which is a function of all oscillation parameters, 
onto any subspace of parameters by local minimization. This approach allows the inclusion 
of multi-parameter-correlations, where external constraints (e.g., on the solar parameters) 
can be imposed, too. Applications of the projection mechanism include the projections onto 
the sin 2 2^13-axis and the sin 2 2#i3-5cp-plane. In addition, all oscillation parameters can be 
kept free to precisely localize degenerate solutions. 

2 The computation of raw event rates 

In this section, we sketch the computation of the event rates, which is one of the core parts 
of the GLoBES software. However, because of the complexity of this issue, we refer to the 
GLoBES manual [14] for details. 
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The differential event rate for each channel is given by 
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where a is the initial flavor of the neutrino, /3 is the final flavor, <& a (E) is the flux of the 
initial flavor at the source, L is the baseline length, N is a normalization factor, and p is 
the matter density. The energies in this formula are given as follows: 

• E is the incident neutrino energy, i.e., the actual energy of the incoming neutrino 
(which is not directly accessible to the experiment) 

• E is the energy of the secondary particle 

• E' is the reconstructed neutrino energy, i.e., the measured neutrino energy as obtained 
from the experiment 

The interaction term is composed of two factors, which are the total cross section a l p{E) 
for the flavor / and the interaction type IT, and the energy distribution of the secondary 
particle k l p(E — E). The detector properties are modeled by the threshold function Tp{E), 
coming from the the limited resolution or the cuts in the analysis, and the energy resolution 
function Vp{E — E') of the secondary particle. 

Since it is a lot of effort to solve this double integral numerically, we split up the two 
integrations. First, we evaluate the integral over E, where the only terms containing E are 
kf(E - E), Tp(E), and V P (E - E'). We define: 

oo 

R, T (E, E') ef(E') = J ' dE T P (E) kf(E - E) V P (E - E') . (2) 

o 

Thus, Rp T (E,E') describes the energy response of the detector, i.e., a neutrino with a 
(true) energy E is reconstructed with an energy between E' and E' + dE' with a probability 
Rf(E, E')dE'. The function R(E, E') is also often called " energy resolution function" . Ac- 
tually, its internal representation in the software is a smearing matrix. The function e l p(E') 
will be refered to as "post-smearing efficiencies", since it will allow us to define cuts and 
threshold functions after the smearing is performed, i.e., as function of E' . In addition, 
GLoBES uses "pre-smearing efficiencies", which are evaluated before the smearing is per- 
formed, i.e., as function of E. Similarly, constant 1 background rates can be added to the 

1 With respect to the oscillation parameters, not the energy. 
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Figure 2: The most important components of AEDL: Channels, rules, and experiments. 

event rates before or after the energy smearing, which are referee! to as "pre-smearing back- 
grounds" and "post-smearing backgrounds". These types of efficiencies and backgrounds 
allow a very accurate modeling of many factors, such as atmospheric or geoneutrino back- 
grounds, energy cuts, or energy threshold functions. 

Eventually, we can write down the number of events per bin % and channel c as 

r E l+ A El /2 d IT 

< = / dE' -JL{Ef) (3) 

where AEi is the bin size of the zth energy bin. This means that one has to solve the integral 

OO 
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which gives the raw event rates of the channel c in the zth bin. Note that the events are 
binned according to their reconstructed energy. 

Core part of the event rate computation is the energy smearing algorithm to evaluate Eq. (0J), 
where either a particular energy resolution function can be chosen for automated energy 
smearing, or the discretized smearing matrix R C (E, E') = can be given manually. In 
addition, it is possible to define a low-pass filter to avoid aliasing effects for neutrino oscil- 
lations faster than the sampling width. 



3 Definition of experiments with AEDL 

(Abstract Experiment Definition Language) 

In order to define experiments, GLoBES uses AEDL ("Abstract Experiment Definition Lan- 
guage"). An experiment normally corresponds to an AEDL file, which is a human readable 
text file written in AEDL syntax. 

The key concept of AEDL is to regard a detector as an abstract system which maps the true 
properties of a neutrino onto the reconstructed properties of the neutrino. The latter are 
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subsequently used in the fit of the oscillation parameters. Within GLoBES, only energy and 
flavor are observables, which is sufficient for long-baseline and reactor experiments. In other 
experiments, such as atmospheric neutrino experiments, more observables (e.g., direction) 
may have to be considered. 

The main components of any AEDL experiment definition are "channels", "rules", and "ex- 
periments" (cf., Fig. |2J). A channel corresponds to a neutrino oscillation channel including 
flux, cross section (for one specific interaction type), energy resolution function, initial and 
final neutrino flavors, their polarity (neutrinos or antineutrinos) , and efficiencies. The chan- 
nel raw event rates are computed according to Eq. (j3J) . Each channel leads to the raw event 
rates for a specific interaction type. In AEDL, many of the channel components have to be 
defined or loaded from files before, such as fluxes, cross sections, or the energy resolution 
function. For the different options, see the GLoBES manual [14]. 

For a "rule" , the raw event rates of one or more signal channels and one or more background 
channels are added. The splitting in signal and background is arbitrary, but all of the signal 
and background components are defined to have a common systematics. The event rates of 
all signal and background components add up to the total event rate of the rule, which leads 
to a (Ax 2 ) r - The signal or background within each rule allows the specification of signal and 
background normalization errors and energy tilt or calibration errors. These systematical 
errors are evaluated with the "pull method" such as in Refs. [7,15]. In addition to these 
systematics errors, an overall evaluation strategy is assigned to each rule, which specifies 
the type of systematics (tilt or calibration error), and the use of spectral information or 
total event rates. 

Finally, one ore more rules add up to an experiment, where the total A% 2 is obtained as 
the sum of the (A% 2 ) r of all rules. This approach allows the definition of appearance and 
disappearance channels, neutrino and antineutrino running, or interaction types with dif- 
ferent systematics (spectral information versus counting rate) within one experiment. The 
GLoBES user interface allows the simulation of one or more experiment simultaneously, 
which means that one could also use different experiments for different oscillation channels. 
However, there is still one component missing in the experiment definition, which is the mat- 
ter density profile. For an experiment, an arbitrary matter density profile can be specified, 
which is evaluated with the evolution operator method [16]. In addition, the matter density 
profile is multiplied by a scaling factor, which is treated as an independent parameter with 
a relative (matter density) uncertainty. With this approach, one can take into account that 
the matter density profile along a specific baseline is only known to about 3% — 5%. Thus, 
for one particular baseline, all rules should be defined within one experiment. 

4 Analysis of experiments: The C user's interface 

In order to use GLoBES, a C-library is linked with the application software. This library 
provides the user's interface functions. It allows to load AEDL files, initialize the experi- 
ments, and have access to various Ax 2 -functions including any combination of systematics 
and correlations. In addition, it provides low-level access to oscillation probabilities and 
event rates, and allows the readout and modification of many experiment parameters at run 
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Figure 3: Left plot: Example for a correlation between sin 2 29\^ and 5cp (for 1 d.o.f. only, un-shown 
oscillation parameters fixed). Right plot: The % 2 -value of the projection onto the sin 2 2#i3-axis as function 
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the position of these minima in the left-hand plot. The arrows mark the obtained fit ranges for sin 2 20\^ at 
the 3cr confidence level (1 d.o.f.), i.e., the precision of sin 2 2#i3. 



time. 

The most sophisticated feature of the user's interface is the possibility to include the full 
multi-parameter correlation among all oscillation parameters. Together with the matter 
densities of the iV ex p loaded experiments (iV exp > 1), the neutrino oscillation parameter 
space in GLoBES has 6 + N exp dimensions. In the x 2_a PP roacri , the allowed region of k 
(1 < k < 6 + iV ex p) fit parameters (at the chosen confidence level) is obtained as the 
projection of the 6 + iV ex p- dimensional fit manifold onto the fc-dimensional hyperplane. 
For example, for the precision of the parameter sin 2 2# 13 only, the fit manifold needs to 
be projected onto the sin 2 26 l 13 -axis. One can easily imagine that the computation of the 
full multi-parameter correlation is very expensive with a grid-based method. However, the 
topology of the neutrino oscillation parameter space is well investigated and rather smooth. 
Therefore, it is feasible to use local minimization for each solution in the parameter space. 
For this process, a 6 + iV ex p- dimensional local minimizer is started at the "guessed" position 
for the solution (such as from analytical knowledge) and then runs into a local minimum. 
This exercise has to be repeated till all local minima, i. e. degenerate solutions, are found. In 
Fig. El this minimization is illustrated at the simple example of a two-parameter correlation 
at a neutrino factory. 

As we have indicated above, the matter density of each experiment is treated as an addi- 
tional parameter, which is known from external measurement with a certain precision. In 
principle, the GLoBES user's interface allows such externally imposed precisions for all os- 
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cillation parameters. This feature turns out to be especially useful for input from the solar 
neutrino experiments, since they provide better precisions for the solar neutrino oscillation 
parameters Am^ and 612 than long-baseline experiments themselves. 

5 Complexity and computational issues 

The GLoBES software only uses polynomial algorithms, which are, however, suffering from 
the large number of dimensions. The Ax 2 computation basically consists of two steps: The 
systematics (Ax 2 ) sys computation with glbChiSys, and the projection algorithm which uses 
glbChiSys many times. The overall computation time is obtained as the product of these 
two steps. 

The run time of the systematics (Ax 2 ) sys is 

tsys O (-^Sampling ^ -^Bins -^Sampling X -/VL ayers ) , (5) 

since the oscillation probabilities for all channels and energies are only computed once and 
stored in a list. Here iVgampiing is the number of sampling points for the evaluation of the 
right integral in Eq. (jlj), -/Veins is the number of energy bins, and -/Vipers is the number of 
matter density layers. Thus, whenever the number of matter density steps is small, the 
product of sampling points and bins dominates. In practice, the iVs amp i ing x iVei ns energy 
smearing matrix is already computed at the experiment initialization in order to save run 
time. In addition, since it contains a lot of zeros, GLoBES uses a special format to avoid 
looping over the zeros. 

The run time for the projected A% 2 would, for a grid based method, be 

tgj = O (n 6+ *-) x t sys , (6) 

where n is the number of grid points to be evaluated for each oscillation parameter, and 
-/V exp is the number of experiments (since the matter densities of all experiments appear as 
parameters). Though it is, in principle, possible to use such a method with GLoBES, it is 
impracticable in most cases, since it would involve at least billions of steps for the complete 
multi-parameter correlation. The local 6+-/V exp -dimensional minimization reduces this effort 
to 

t MiBimi Z er = ( mQ) x ^ x ^ (7) 

steps, where -/Voeg is the number of (disconnected) local minima, i.e., the number of degen- 
eracies (typically O(10)). 

Eventually, GLoBES needs about 10 to 15 seconds to compute a projected A% 2 on a modern 
Pentium machine. Therefore, more sophisticated applications, such as a two dimensional 
plot for CP violation as function of the simulated values of sin 2 26 Vi and <5cp can be obtained 
in several hours. 

6 Use of package 

The GLoBES package [17] is a tar-ball which consists of the source code for the C-library 
providing the user's interface, the source code for a program to test AEDL files, example 
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C-files illustrating the use of the library, experiment prototypes in AEDL with supporting 
flux and cross section files, and the usual supporting files for installation and compilation. 
Currently, GLoBES is only provided for GNU/Linux. In addition, an extensive manual 
covering the user's interface and AEDL [14] can be obtained from the web-site. GLoBES is 
free software and as such licensed under the GNU Public License. 

The installation of GLoBES is highly automated by the use of the Autotools family and 
follows the usual procedures. By default a shared library libglobes.so is installed, but a 
static version is available, too. In addition, an executable named globes is installed. It is 
another possibility to install GLoBES without root privilege into a user's home directory. 
During the automatic configuration process also an example Makefile is produced which can 
serve as a template for compiling and linking own applications against libglobes . so. This 
Makefile can be used to directly compile and link the examples from the manual. 

AEDL files, such as the experiment prototypes, can be edited with any text editor, in order to 
be loaded by the users's interface later. In addition, the globes binary allows to develop and 
test experiment descriptions. For example, event rate information can be quickly provided 
to adjust the neutrino flux normalization. For further information, we refer at this place to 
the GLoBES manual [14]. 

7 Summary and conclusions 

In summary, the GLoBES software package provides powerful tools to do a full three-flavor 
analysis of future long-baseline and reactor experiments including systematics, correlations, 
and degeneracies. In addition, the abstract experiment definition language allows to define 
experiments at an abstract level in a highly accurate fashion. Some of the major strengths 
of GLoBES are the ability to quickly define new experiments, the potential to take into 
account the full multi-parameter correlation, the possibilities to include external input and 
matter density uncertainties and the ability to analyze several experiments simultaneously. 

We conclude that the GLoBES software has two major target groups: Experimentalists, 
who want to quickly evaluate the physics potential of their setups, and theorists, who want 
to test new ideas or strategies with pre-defined experiments. Especially, the separation 
between experiment descriptions as simple text files (together with supporting files) and 
the application software should allow an efficient and lively interaction between those two 
target groups. 
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